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I/-S ■ Abstract 

OO , Vortex solutions in 17(1) x f/(l) Chern-Simons theory coupled to a pair 

^3 ■ of hard-core bosons representing two layers of electrons are analysed. It 

l^ ' is shown that there is such a range of parameters (a/3 < 7^) in which 

finite charge (1,0) or (0,1) vortices can not exist. Instead the minimal 
flux quanta are (1, 1) hybrids. Their statistical dominance is shown to 
have an influence on Haldane-Halderin hierarchy of states. The form of 
Coulomb potential and its influence on the range of physically acceptable 
states is discussed. 
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^ ■ 1 Introduction 

O 

O ' The two-dimensional electron system can have a rich structure. One of the 

most significant discoveries is the fractional Hall effect in strong magnetic fields 
- [nl . The standard theoretical approach is through variational analysis with trial 

^S ' functions M. There are also ficld-thcoretical studies of the subject based on 

^ . generelisations of Ginzburg-Landau model with a use of effective Chern-Simons 

theory [BJ. More microscopical treatment of the phenomenon was introduced 
by Ezawa and Iwazaki B . Recently there was discovered the even-denominator 
filling factor in double- layer electron system H . The system was proposed to be 
described by U{1) x U{1) Chern-Simons theory coupled to a pair of hard-core 
bosons representing two layers of electrons . Nonperturbative excitations are 
topological vortices with some winding numbers on both of the layers {p, q) . If 
(1, 0) and (0, 1) vortices condense the hierarchy of states was shown to begin 
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with characteristic iq instead of 2g as in a single-layer system (q is an integer). 
By carefuU analysis we show that there is such a range of statistical parameters 
in which minimal vortex excitations with finite charge are (1,1) vortices living 
on both of the layers. In this respect this work is a generalisation of our previous 
paper |^ . The physical consequence of this fact is that the denominator is once 
again 2q. 

Another question is the form of Coulomb potential. To have finite Coulomb 
energy of the condensates one needs regularising background. If the background 
charge density is just opposite to the condensate density on a given layer the 
whole spectrum of states (k,l,m) is allowed. If instead the background is uni- 
form then only states with k — I survive switching on the Coulomb interaction. 

2 The model, condensates and vortex solutions 

The model we consider in this paper was introduced by Ezawa and Iwazaki . 
It consists of two species of polarised electrons represented by hard-core bosons 
coupled to U{1) X C/(l) Chern-Simons field the effect of which is via singular 
gauge transformation equivalent to imposing the Fermi-Dirac statistics. 

_ ^ M>>A (l)g (1) _ J_ A^-^A (2)^ (2) 
^-^^ "m ''A ^Y t^ A 

The coefficients X, Y, Z are defined in an entangled way by another set of pa- 
rameters 

X= "^-^' , y^ «/^-7^ Z^ l'-^P (3) 

p a 7 

These parameters take values a = nk, (3 = ^1,^ — ttto with k,l,m odd integers 
which provide electrons ■0: 4' with desired statistics. The potential term has 
the form characteristic for bosonic end perturbation theory with respect to the 
parameters a,/3,7, see Q] and e.g.[g| 

V=^[a\^\'+2^\^\'\^\'+P\<Pn , (4) 

This potential is a field-theoretical remnant of the two-particle delta-type po- 
tential g. Another potential term is a nonlocal Coulomb interaction 

AV'^ = AV^^ + AVg + AVg . (5) 

The particular terms are given by 

An'i = iJ / ^'^ d'y t^x^) - Pi] ^^ [^*VXy) - pI] , 
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where e is a dielectric constant and d is the interlayer distance. We have taken 
into account uniform background the charge density of which is opposite to the 
medium charge density of electrons. The Lagrangian in this form does make 
sense only if a/3 — 7^ 7^ - the case we refer to as nondegenerate, but the model 
itself can be made more general. Couplings between electrons and gauge fields 
are established by 

i?^V = S^V' - «<V + «eA^"V , D^^^d^4,^ia''^U + ieA1-'^- (7) 

In what follows we will consider only external magnetic field and put Aq^^ = 0. 
Variation with respect to Lagrange multipliers % , 7=1,2 leads to constraints 

Bi = 2api + 27P2 , 

B2 = 27pi + 2/3p2 , (8) 

where the electronic densities are pi — ip*'ip , p2 — (j)*4'. The Hamiltonian of 
the model is 

For the following we will switch off the Coulomb interaction and later on regard 
it as a correction. In a uniform external magnetic field Bogomol'nyi decompo- 
sition Igl leads to 

/I 1 ^ ryext 

'^''^ [ 2M ' ^^' + *^'^^' '' +2M ' ^^' + '^'^'^ '' ^^2Nr^P' + P''^ 

+V - ^{PiBi + P2B2) - V X Ji - V X J2 ] . (10) 

If we neglect boundary terms and take into account Gauss' laws (|8|) the second 
line will vanish identically. The term pi + p2 in the first line can be integrated 
out if we impose an extra costraint of definite number of particles iV 0, 0| . Thus 
we are left with the Hamiltonian 

H^\..N + jd^.[±-^\D,^\^^±-^\D,^n^ (11) 

where u>c = '^\j is a cyclotron frequency. Minimal energy static solutions 
satisfy self-dual equations 

D+t/- EE a+^ - ia'^^ V + ie^T*^ = ' -D+0 = , (12) 



in addition to Gauss' laws (g). Solutions to these equations are also solutions 
of full Euler-Lagrange equations of the model only if Lagrange multiplyers are 
equal to 

% = -api - 7P2 , flo = -7Pi - PP2 ■ (13) 

For a solution with constant densities pi , p2 the phases of scalar fields can be 
put constant and Eqs.(|l^) imply that 

Bi = eB"^^ , B2 = eS'^^* . (14) 

These and the constraints (g) lead to the following values of uniform condensates 
in the nondegenerate case 
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For positive value of i?^^*, which we take for definiteness, the nonzero conden- 
sates exist if and only if a — 7,/? — 7 and a/3 — 7^ are simultaneously positive or 
negative. The total electronic filling factor in this case is 



i^ 



(0) _ 27r(pi +P2) _ TT{a + l3-2-/) _ k + I - 2m 
ei?'^^* aP — 7^ kl — m? 
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Thus for e.g. (k,l,m) = (3,3,1) we have the experimentally observed filling 
factor i II . 

In the degenerate case there is a continuous family of constant solutions 
characterised by a constraint 

p = p^ + p^ = —— , (17) 

27 

only if in addition a = /3 = 7. The filling factor is equal to i^*-"-* — — with m 
being an odd integer. A free parameter in the solution signals existence of a 
zero mode. 

Thanks to the form of Coulomb potential all the uniform solutions (g|) have 
zero electrostatic energy. On the other hand contrary to the assumption in (|6|) 
one could think it to be more natural that the positive charge density of the 
background is constant everywhere and in particular the same for the two layers 
yO° = i(p° + P2)- I'^ such a case the Coulomb potentials should be taken as 

2e J \x-y\ 

AVg = - fd'x dPy [r^{x) ~ p°] . \^ ^^ [0>(y) - p'] , 
e J ^{x - y)2 + (P 

^y?i = f /^'^ d^y W^{x) - p^\ y^— [0>(y) - p^\ . (18) 

2e J I a; ~ 2/ I 



With this background any constant solution such that pi ^ p2 would have 
infinite energy. Thus the only physical states would be those with a — ^ 
or in other words one could observe only states of the type {k,k,m), which 
characterise by filling factors z^'"' = -, where s is any integer (even or odd). In 
the degenerate case there would be no longer any free parameter but rather the 
separate charge densities would be locked at pi = P2 = ^ — • Thus the zero 
mode can be expected to become massive thanks to Coulomb energy as it was 
observed experimentally g . 

Now let us take into account a possibility of nonuniform vortex solutions. 
Let us write the scalar fields in the form Tp = y^expiwi, (/) = Y^expia;2- 
The self-dual equations ( |l2|) combined together with Gauss' laws (||) lead to the 
following modified Toda equations 

iv^lnpa = 27P1 + 2/3p2 - eS""* + em„9m5„W2 , (19) 

The phases luj are single- valued except of finite sets of singular points. Outside 
of these points one can neglect the last terms in above equations. It is easy to 
see that if we impose a condition p2 = §E^Pi inspired by the uniform solutions 
5|) we will obtain an equation 

iv^lnpi = ^^^^Pi-ei.-, (20) 

2 p-7 

which is known to possess multivortex solutions. Thus in any case there are 
static multivortex solutions with vortices with winding numbers (1,1) located 
at arbitrary points of the plane. We will call them hybrid anyons. Now the 
question is whether we can also expect solutions with separate vortices of the 
types (1, 0) and (0, 1). Clearly if we put artificially 7 = then eqs.(pj|) decouple 
from each other and become similar to (pG|), so in this limit there exist such 
isolated vortices. The opposite limit oi a — /3 = and 7^0 characterises 
by existence of only hybrid (1,1) vortices as physical excitations as was shown 
in |7|. In a more general case we have to analyse carefully asymptotices of 



such potentially interesting solutions. We substitute to Eqs.(|l9[) following forms 
Pi — p^{i + fi) and linearise them with respect to //'s 



ial3 - 7')V'/i = 2a{l3 - j)fi + 2j{a ~ 7)/; 



2 



{aP - 7')VV2 = 27(/3 - 7)/i + 2/3(a - 7)/2 , (21) 

Of course there is a solution with /i = /2 which corresponds to the exact 
solution with only hybrid vortices. To look for a more general solution let us 
take f2 = fi+ u: 

ap — 7^ 



27(a-7 ). 
aP — 7^ 



VVi - 2/1 ^ ^fT^" > (22) 



The first equation is an eigenvalue problem which is solved by an asyniptotics 
of Bessel or modified Bessel function dependent on whether the eigenvalue on 
R.H.S. is negative or positive respectively. Its sign is the same as the sign of its 
denominator a(3 — 7^ . If this determinant is negative then both /i and /2 have 
an admixture of slowly vanishing solution 

^1 = ^ - — rP^" ' -^2 = /+ f:^^ u , (23) 

a + p — 27 a + p — 27 

where / is a general exponentially vanishing solution of uniform equation V"/ — 
2/ = 0. The charge density fluctuations due to u are 

A. ^ pext («"7)(/?-7) r r /o/l^ 

dpi = -T^eB —- — — — — -—u , 5p2 = -dpi . 24 

2 (ap — 7^)(a + p — 27) 

A total electric charge of such a vortex does not have any contribution from these 
fluctuations since they exactly cancel each other. Nevertheless separate charges 
on each of the layers are divergent and not well localised. This signals that if 
we take into account Coulomb interaction such solutions even if substantially 
distorted can have much higher energy then (1, 1) vortices. 

For a positive determinant u vanishes exponentially and there are no prob- 
lems with charge normalisation. Thus we can see that the line in the parameter 
space a(3 = 7^ is a border between an area where as physical excitations there 
are only hybrid (1,1) vortices (a/3 < 7^) and that in which there is a possibility 
of flnite charge and energy separate (1,0) and (0,1) vortices (a/3 > 7^). In 
more physical terms if the interaction 7 between the two layers is large enough 
as compared to a and (3 vortex can not form only on one of the layers. Only 
the double-layer (1, 1) composites are allowed. 

Similar analysis can be performed at the bifurcation point a = /3 = 7. Eqs. 
( p2[ ) are replaced by 

V^w^O , 

VVi ~h^\u ■ (25) 

A general solution is given by /i^2 = / + o^o- uq is a constant which is a 
parameter connected with the zero mode - it only changes asymptotic values 
keeping p\ + p2 flxed. This symmetry is broken and the parameter is locked at 
uo = if we take the Coulomb interaction of the form (|l8| ) . / is an exponentially 
vanishing asymptotics of modifled Bessel function. Thus if a = /3 = 7 and both 
condensates are nonvanishing the smallest quasiparticle excitations are hybrid 
(1, 1) vortices. This perturbative result is in agreement with conclusion in pO| , 
Eq.(2.19)] if we take into account that in their paper boundary conditions are 



imposed in such a way that one of the condensates vanishes. Our boundary 
conditions can be obtained by global SU(2) transformation. 

Here we have just analysed properties of eventual single isolated rotationally 
symmetric (1,0) or (0, 1) vortex solution. Now we will try to directly compare 
properties of (1, 1) vortex with a pair of (1, 0) and (0, 1) vortices. 

3 Splitting of hybrid vortex into separate (1,0) 
and (0, 1) vortices 

Keeping in mind the whole area aP < 7^ we will restrict in this section for a 
sake of simplicity to the case a ^ f3 — and 7 > 0. Let the unperturbed hybrid 
vortex be 

1/; = = FWe*" . (26) 

Now we take perturbations of the phases of the scalar field to be ai,a2 and 
those of the moduli: Fhi , Fh2- 

+ 50 = i^(r)[l + /i2(r,0)]e'"'+*"=('-'^) . (27) 

Linearisation of the self-dual equations (03) with respect to perturbations of the 
scalar fields and those of gauge potentials Cj, , / = 1, 2, yields 

Ca — drhj deaj , 

r 

4''' ^ --dghi ~ drai (28) 

r 

for / = 1,2. Once the perturbations of the Higgs field are known, c). can be 
calculated from the above equations. To have a unique solution we have to fix 
the gauge 

dkci'^ = dj/^ + ic(^) + ^declf^ = , 7 = 1,2. (29) 

Upon substitution of (ESI) to this gauge condition we will obtain 

V^ai = , V^as = . (30) 

Similar substitution of ( p8| ) to Gauss' laws (^) will lead to 

V^hi=A-fF^{r)h2 , V^h2^-i-fF^{r)hi . (31) 

We would like to describe decay of the hybrid vortex in the direction of x-axis 
on two separate vortices. That is the reason for symmetry requirements on hi 

hi{x,y) ^h2{-x,y) , hj{x,y) = hi{x,-y) . (32) 



They restrict Fourier transforms with respect to 6 to foUowing forms 

hi (r, e)^Ha{r) + Hi (r) cos 6i + (5/ii , 

/i2 (r, 61) = Hq (r) - Hi (r) cos 9 + Sh2 , (33) 

where 6hi are tails of Fourier seria beginning at terms cos 20 and sin 20. Sub- 
stitution of these seria to Eqs.(|3l|) leads to equations 

{fj + ~)Ho^'i^Ho , 
dr'^ r dr 

(;^ + ~- >.—•-,//,, (34) 

Thus Hq is an exponentially vanishing function while an asymptotics of Hi is 
equal to zero or proportional to the asymptotics of Bessel function and slowly 
vanishing. But why should Hi be nonzero? 

Its asymptotics near r = can be read from Eq.(B4), namely Hi{r) ^ ^. 
With such an asymptotics the scalar fields for small r = and small ^ tend to 

^^re*e_^ , 0-re*''+^ , (35) 

and zeros are shifted to ^ and — ^ respectively. Thus if (1,0) and (0, 1) vortices 
are already not exactly on top of each other then ^ ^ and also _ff(i) has to be 
nonzero. 

Putting all the above together let us calculate changes of electron numbers 



SQi = d'x[ F\r)il + hi)' ~ F'{r) ] . (36) 

With the forms ( |33| ) the above integral can be finally rewritten as 

SQi = f d^x F^{r)[2Ha+Hl] + j d^x F^{r)[H^ + {Shif] . (37) 

The first integral is finite. To show that the second one is infinite it is enough 
to calculate its first term 

I d^xF^{r)Hl{r)^ constR , (38) 

where i? is a radial cut-off. Thus we can see that JQ/'s are badly divergent. 
This divergence is not an artifact of linearisations we have done since it depends 
only on the asymptotics of -ff(i) for large r where the approximations are valid. 
In other words the lowest energy hybrid (1,1) vortex can not decay into sep- 
arate (1, 0) and (0, 1) vortices because such a decay would violate conservation 
of the numbers of particles. The decay could be possible only from an excited 
state of hybrid vortex. A pair of separated (1,0) and (0, 1) vortices must be 
higher energy configuration than a single hybrid vortex. Once again we are lead 
to conclusion that hybrids will dominate statistics of the system. 



4 Condensation of vortices and hierarchy of 
states 

Let us consider the nondegenerate case a(3 — 7^ 7^ 0. We perform a singular 
gauge transformation 

uji ^ uji + ai , a^^) -^ a'^J^ + d^aj , (39) 

where positions of vortices are encoded in singular shifts of phases 

a7-^e(z-z,,) , (40) 

Si 

on the fields in the Lagrangian 

L = t^/Dn^/j - —Dk^*Dk^lj + ir ^00 - —Dk<j)*Dk<p ~ V{il;, 0) 

- ^e^'''4'dJ^^-^e^'^>^a^^dJ^^-V , (41) 

Under such a transformation which is intented as an introduction of new vortex 
degrees of freedom the Lagrangian changes by 

^-'-'vortex — ^~[ ■^ 1 -rp ' y ) 

1 d^aiK'^ d^aiK'^ df,aiK!^ + d^,a2K^ 
- -^-i—^- + -^ + z ) ^ (42) 

where vortex currents are 

K^ - ^e'^-'^d^dxaj - J2 is/*'H^ -^sr) , (43) 

si 

with the parametrisation of vortex positions by x^^ = {t,x'^j). With this form 
of the vortex current one can easily rewrite the change in the Lagrangian ( [42[ ) 
as 

lALyortex = — 7r( — \ — 1 j 

-^[^ J2 e(zp, -zsj + ^ J2 0(^P2-^.2)+ ^E.i..,0(^^i-2s2)](44) 

Pl>Sl P2>S2 

The charges of vortices with respect to the C-S fields can be read from constraint 
equations 

Api = ABi AB2 , 

AP2 = -^Ai?i - ^AS2 , (45) 



where by A we denote deviations with respect to uniform background. Thus 
upon integrations and use of quantisation conditions we find that charges of 
vortices with winding numbers (p, q), if they exist, are given by 

(3p~'jq 

qi[p,q\ =7r— , 

ap — 7^ 

aq--fp 

92P,g =Tr—5 2 ' ('^6) 

ap — "f^ 

while their electric charges are 

eQiP, q\ = ~~e{qi + q2) = -ctt ^ . (47) 

ixp 7 

Let us concentrate on the case a/? — 7^ < in which dominating vortices 
are those with winding numbers (1,1) - the hybrid anyons. The vortex currents 
can be identified K^ — K^ — A'*', the charges gi = tt J^T^ ^ and (72 = tt ^^^ ^ 
are positively definite. The vortex part of the Lagrangean can be rewritten as 

p<q 



where the statistical parameter k is defined modulo integer multiplicity of 2tt 
K = -n ( 2~)~27rg . (49) 



Upon variation of the effective Lagrangian with respect to Og we get a new set 
of constraint equations 

Pi = ^ 2X^^ ~ 27"^^ ^ qiPvortex , 
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P2 = ^^^1 ~ oy -^2 + q2Pvortex ■ (50) 

With the identification K'^ = Pvortex these equations are to be understood in 
terms of mean values. For large enough separations of vortices the Lagrangian 
(ksf) can be supplemented by kinetic term similarly as in M 



AiJ™We. = ^E[Pp-9l4'^-924'^+Cfc]' . (51) 

" P 

The auxiliary vector field is given by Ck{z) = ^ J2p dkQ{z—Zp) or by a constraint 



It is worthwhile to mention that the form of the Hamiltonian (51) is justified 



only if separations of vortices are large as compared to sizes of their cores. In the 
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derivarion vortices were represented by infinitely thin delta-like distributions of 
magnetic flux and charge. Short range interactions of hybrid (1, 1) anyons were 
analysed in H in the case of only mutual interactions. The conclusions can be 
easily generalised to this case. The main result is that statistical interactions 
mediated in (^ by the fields Cfe are effectively switched off at short distances. 
Instead there are short range charge-flux interactions as if vortices were finite- 
width charged selenoids. They can lead to magnetic trapping and quantum 
mechanically to scattering resonances characterised by a set of Landau levels. 
The Hamiltonian (|5^) leads to the following field-theoretical description 

Hvortex = :rjT- \ {Ok " ^^lOfc " ^20,^ + '>''^k)u ^ , (53) 

The self-dual equation for vortex field D^u = in a uniform state implies that 
qiBi + 52^2 + SmndmCn — 0. This and the relations Bi = B2 — eB^xt lead to 
the filling factor for vortices 

^'^Pvortex ^ fKA\ 

^vortex ^ TTd — ~ ' \ / 

Q-Dext K 

From equations (|5C|) we get the electronic filling factor 



(1) = 



1 . 1 



i^' ' = 



- ' ^=7757 • (55) 



e+^ ' ^ -''' 



This formula is the same as in the single layer quantised Hall effect H . As such 
it has to be compared with that for the double-layer effect when there are single 
vortices of the types (1, 0) and (0, 1). Such a formula was derived in the article 
by Ezawa and Iwazaki p]. 



ly 



(1) 



1 1 



'' ' 4(p+r) ^ ^ 4g 



(56) 



where p, r are integers and can be replaced by q — p + r. The characteristic 
difference between formulas (pq) and ( jSq ) is the multiplyer 2 instead of 4 when 
there is dominance by hybrid vortices (1,1)- 

5 Conclusions 

There are two characteristic areas in the parameter space. If a/S > 7^ vortices 
with winding numbers (1,0) and (0,1) dominate while for af3 < 7^ quite as 
well as for the bifurcation point a = /3 = 7 the most important type are (1,1) 
hybrids. This dominance by one type or the other can be observed thanks to 
characteristic changes in Haldane-Haldcrin hierarchy of states. 

The hybrid (1,1) vortices can be thought of as composites of separate (1,0) 
and (0, 1) vortices. There is always present charge-flux interaction between them 
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being a remnant of long range statistical Aharonov-Bohm-type interaction. This 
interaction acts also between whole hybrid vortices and can lead to magnetic 
trapping and at the quantum level to scattering resonances characterised by 
Landau levels. If a(3 < 7^ then (1, 0) and (0, 1) vortices are additionally at- 
tracted by potential force to form a hybrid vortex. 

The distribution of positively charged background may be important for the 
range of physical states. If the background is the same for the two layers the 
Coulomb interaction will allow only states of the type (fc, k, m) with two charge 
densities equal one another. 
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